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Fig. 1. (Leit) Dictyostelium discoideum. The lower edge, -with the 
smooth region, is the front side of the cell. Its length is of the 
order of 10 ^m. (Right) Experimental trajectories, x{t) — x{0) 
(solid line) and y{t) — j/(0) (dotted line), are plotted as functions 
of time in the case _E = 10 V/cm. 
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Directional cell movement is a fundamental phe- 
nomenon exhibited by many biological processes. It has 
been known that an electric field exists on the surfaces of 
the tissues of organisms and that it acts as a directional 
cue for the type of the cell migration known as electro- 
taxis. Electrotaxis is thought to play important roles in 
various physiological processes including embryogenesis 
and wound healing, and the underlying molecular mech- 
anisms of electrotaxis are now extensively studied. 

In order to elucidate the mechanisms of the electro- 
tactic responses of cells, the cellular slime mold Dic- 
tyostelium discoideum (see Fig. 1 (Left)) is a suitable 
organism to study, because of its high motility and 
strong electrotactic response. With well-established ge- 
netic engineering techniques and advanced microscopic 
techniques, the input-output relationship in the electro- 
tactic response of Dictyostelium cells has been investi- 
gated to elucidate the stochastic processes involved in 
the signaling systems responsible for cell motility and 
their regulations. 

In this paper, we analyze the electrotactic movement 
of Dictyostelium discoideum from the viewpoint of non- 
equilibrium statistical mechanics. Because we can ob- 
serve fluctuating behavior of cellular trajectories, we an- 
alyze the probability distribution of the trajectories with 
the aid of the fluctuation theorem. Recently, the validity 
of the fluctuation theorem was verified in a colloidal sys- 
tem,^) and it has also been applied to granular systems, ''^^ 
turbulent systems^-* and chemical oscillatory waves^^ to 
investigate some of their statistical properties that are 
not yet completely understood. Noting that the fluctua- 
tion theorem is potentially applicable to cellular electro- 
taxis, here we employ it to help us obtain a phenomcno- 
logical model of this biological system. 

System. In this study, Dictyostelium discoideum, Ax2 
cells (wild type) were starved for up to 4 hours, with a 
pulse of 100 nM cAMP applied every 6 min. The cell 
suspension was injected into the chamber for electrotac- 
tic assay, and the cells were allowed to spread over the 
coverslip for 20 min at T = 294 K. Direct current was 
applied to the chamber as illustrated in the paper of Sato 
et al.'^^ 

The cells in the chamber were observed with a micro- 
scope capable of producing differential interference con- 
trast optics. Data acquisition started 5 min after the elec- 
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Fig. 2. (Top) The probability distribution, P{A), in the case 
-E = 10 V/cm with r = 10 sec (filled circles) and t = 15 sec 
(open circles). 4000 samples of A were obtained from 58 inde- 
pendent trajectories. (Bottom) The quantity ln[P(yl)/P(— A)] as 
a function of A in the case r = 10 sec (filled circles) and the case 
r = 15 sec (open circles). The slope of the straight line is equal 
to 2 iim~^. 



trie field was first applied, and the electric field remained 
at a constant strength, E, throughout experiments. Un- 
der these conditions, we considered the system to be in 
a steady state. To analyze the motile activities of the 
cells under the electric field, cell images were processed 
automatically with a time resolution of 5 sec and con- 
verted into binary images by selecting an optimal value 
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of the brightness threshold. In this way, the trajectory 
of the center position of the observed ceUular region, 
{x{t), y{t)), was determined. Because the gradient of the 
electric field is non-zero only along the x direction, we 
particularly investigate the trajectories in the x direction 
for the case that £^ = 10 V/cm, which is sufficiently large 
for the response' saturation of the cell.^^ 

Fluctuation theorem. When the electric field is turned 
on, Dictyostelium cells begin to migrate toward the cath- 
ode (the +x direction). In Fig. 1 (Right), we plot example 
trajectories, x{t) — x{0) and y{t) — y{0), as functions of 
time in the case E = 10 V/cm. The fluctuating behav- 
ior of the cell movement was observed. Then, using these 
fluctuating trajectories and setting A = x{t) — x{0), we 
investigate the probability density for the realization of 
the value A. In Fig. 2 (Top), we plot the probability den- 
sity, P{A), in the case = 10 V/cm with r = 10 sec and 
T = 15 sec. We flnd that the graphs have exponential 
tails. 

In a colloidal system under a non-equilibrium condi- 
tion, an entropy production, S, of the system was mea- 
sured experimentally,'^'' and the probability distribution, 
P(I]), was analyzed. This P(S) also has exponential tails, 
and it is known that the entropy production fluctuation 
theorem ln[P(S)/P(-E)] = S holds in the system. Be- 
cause the forms of the graphs in Fig. 2 (Top) are sim- 
ilar to those observed in the colloidal system, wc plot 
\n[P{A)/P{-A)] as a function of A in Fig. 2 (Bottom). 
In Fig. 2 (Bottom), we obtain the linear slopes in the 
cases T = 10 sec and r = 15 sec, and they differ slightly. 
We thus conclude that the value of the slopes is almost 
independent of r at least to the extent of this rough treat- 
ment. Prom the fitting of the data to a linear function, 
we obtain the value of the slopes 2 /zm"^. 

In order to understand a physical meaning of the value 
2 nm~^, we employ a model simplified enough, noting 
that T is loss than the persistence time of the cell move- 
ment. We consider that a description of the cell move- 
ment is roughly described by the Langevin equation 



dt 

{Emit')) 



2arS(t-t'), 



(1) 



where F is the effective friction coefficient, F is the appar- 
ent driving force, and a is a parameter that characterizes 
the fiuctuations of the cell movement. 

A physical meaning of the value 2 /im~^, which is al- 
most independent of r, is understood by one expression 
of the fluctuation theorem. For the model (1), we can 
derive the relation 



In 



P{A) 



-A. 



(2) 



P(-^) 

The relation (2) means that we can infer the value of F/ a 

by computing P(j4), P(— A) and A from the experimental 
data,^) and that the value of the slope 2 /zm"^ in Fig. 2 
(Bottom) corresponds to F/a. 

Discussion. Because the traction force of a Dic- 
tyostelium cell has been measured experimentally and 



found that it is of the order of 1 nN,^) F is also considered 
to be of the order of 1 nN. Then, using F/a = 2 /xm~^, 
a is estimated to be 0.5 nN^um (= 10^ /cbT*^)- Using the 
relations F {v) ^^^^ = F and {v) ^^^^ = 0.1 /im/sec, F is 
also estimated to be 10 g/sec. Here, v{t) is the velocity 
of a cell defined by v{t) = {x{t + At) - x{t))/At where 
At = 5 sec, and ( )^ denotes the steady state average. 

In this paper, we obtained the rough estimation of the 
parameters F, F and a that characterize the cell move- 
ment under electrotaxis by using the fluctuation theorem 
(2), which is derived for the simple model (1). To confirm 
the validity of the characterization, the above parameter 
values should be compared with those of other organ- 
isms, e.g. Amoeba proteus.^^ Also in future, the presented 
method should be extended to satisfy more precise de- 
scriptions of cell movement, noting that the fluctu- 
ation theorem can be derived for more generic stochastic 
processes. 

In general, theories of non-equilibrium statistical me- 
chanics, such as the fluctuation theorem and the Jarzyn- 
ski equality, have been applied to small non-equilibrium 
systems, such as colloidal systems, RNA systems and 
protein systems. However, as there are few studies in 
which these theories are applied to cell systems, we hope 
that the present study will help lead to exploration in 
this new direction. 
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